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ABSTRACT 

Recently, Cai and Su [Phys. Rev. D 81, 103514 (2010)] found that the sign of interaction Q 
in the dark sector changed in the approximate redshift range of 0.45 < z < 0.9, by using a model- 
independent method to deal with the observational data. In fact, this result raises a remarkable 
problem, since most of the familiar interactions cannot change their signs in the whole cosmic history. 
. Motivated by the work of Cai and Su, we have proposed a new type of interaction in a previous 

work [H. Wei, Nucl. Phys. B 845, 381 (2011)]. The key ingredient is the deceleration parameter q 
in the interaction Q, and hence the interaction Q can change its sign when our universe changes 
from deceleration (q > 0) to acceleration (q < 0). In the present work, we consider the cosmological 
constraints on this new type of sign-changeable interactions, by using the latest observational data. 
We find that the cosmological constraints on the model parameters are fairly tight. In particular, 
the key parameter /3 can be constrained to a narrow range. 
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I. INTRODUCTION 



In the dark energy cosmology the well-known cosmological coincidence problem has an important 
position. This problem asks: why are we living in an epoch in which the densities of dark energy and 
matter are comparable? Since their densities scale differently with the expansion of our universe, there 
should be some fine-tunings. To alleviate this coincidence problem, it is natural to consider the possible 
interaction between dark energy and dark matter in the literature (see e.g. I n fact, since the 

natures of both dark energy and dark matter are still unknown, there is no physical argument to exclude 
the possible interaction between them. On the contrary, some observational evidences of the dark sector 
interaction have been found recently [l3L H3 |. In the literature, it is usual to assume that dark energy 
and dark matter interact through a coupling term Q, according to 

p m + 3Hp m = Q , (1) 
pde + ZH(p de + p de ) = -Q , (2) 

where p m and p de are the densities of dark matter and dark energy (we assume that the baryon component 
can be ignored); p de is the pressure of dark energy; H = a/a is the Hubble parameter; a = (1 + z)^ 1 
is the scale factor (we have set ao = 1; the subscript "0" indicates the present value of corresponding 
quantity; z is the redshift); a dot denotes the derivative with respect to cosmic time t. Note that Eqs. (fTJ) 
and ([2]) preserve the total energy conservation equation p tot + 3H(p to t +Ptot) — 0, where ptot = Pm + Pde- 
Since there is no natural guidance from fundamental physics on the interaction Q, one can only discuss 
it to a phenomenological level. This is the realistic status of the interacting dark en ergy m odels so far. 

The most familiar interactions extensively considered in the literature (see e.g. 0t12|) include Q = 
3aHp m , Q — 3f3Hp to ti and Q = 3r]Hp de . It is easy to see that these interactions are always positive or 
negative and hence cannot give the possibility to change their signs in the whole cosmic history. However, 
recently Cai and Su [l5[ found that the sign of interaction Q changed in the approximate redshift range 
of 0.45 < z < 0.9, by using a model-independent method to deal with the observational data. Obviously, 
this result raises a remarkable problem. Motivated by the work of Cai and Su, we have proposed a new 
type of interaction in a previous work , which is given by 

Q = q(ap + 3/3Hp) , (3) 

where a and j3 are both dimensionless constants; the energy density p could be p m , ptot and p de for 
examples; the deceleration parameter 

" H 



^ ~ aH 2 H 2 

Obviously, this new type of interaction Q can change its sign when our universe changes from deceleration 
(q > 0) to acceleration (q < 0). In fact, the deceleration parameter q in Q is the key ingredient of this 
new interaction, which makes our proposal different from the previous ones considered in the literature. 
Note that the term ap in Q is introduced from the dimensional point of view (we refer to fl(| for details). 
One can remove this term by setting a — 0, and then Q becomes simply Q = 3(3qHp (in fact this is the 
very case which will be considered in the followings). 

Since the appearance of the deceleration parameter q in the interaction Q looks speculative to some 
extent, we would like to say some words before going further. Firstly, as is well known, in the literature 
there is no natural guidance from fundamental physics on the interaction Q, one can only discuss it 
to a phenomenological level. In this sense, the other familiar interactions extensively considered in the 
literature have no better origin from the fundamental physics than the one proposed in Eq. ([3]). Secondly, 
we note that q = — 1 — H/H 2 from Eq. (0]) and H 2 oc ptot from the Friedmann equation. Thus, one can 
regard the deceleration parameter q = f(ptot> Ptot) as a function of the total energy density ptot = Pm+Pde 
and its derivative. In this sense, the interaction Q — q{ap + 3f3Hp) = f(p,p) is not so unusual, since it 
is reasonable to image that Q depends on the energy densities of dark energy and matter. Finally, while 
the familiar interactions extensively considered in the literature (such as Q = 3aHp mi Q = 3(3Hp to t, and 
Q = 3i]Hp de ) cannot give the possibility to change their signs in the whole cosmic history, our proposal 
in Eq. ^ provides a possible way out. So, we consider that it deserves further investigation. 
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In [l6|, we have studied the cosmological evolution of quintessence and phantom with this new type 
of sign-changeable interactions, and found some interesting results. In the present work, we would like 
to consider the cosmological constraints on this new type of sign-changeable interactions, by using the 
latest observational data. For simplicity, in this work, we restrict ourselves to the decaying A model (see 
e -g- [13 an d references therein), namely, the role of dark energy is played by the decaying vacuum energy. 
In this case, Eq. @ becomes 

PA = —Q ■ (5) 
The Friedmann and Raychaudhuri equations are given by 

K 2 K 2 

H 2 = -jPtot = y(PA + Pm) , (6) 

K 2 K 2 

H = — -(ptot+Ptot) = — 2 Pm ' ^ 

where k 2 = 8ttG. Notice that we consider a flat Friedmann- Robertson- Walker (FRW) universe throughout 
this work. In Sec. [TTl we briefly introduce the latest observational data which will be used in this work. 
In Sec. IIII1 we consider the cosmological constraints on three particular sign-changeable interactions, i.e., 

Q = q(ap m + 3(3Hp m ) , (8) 

Q = q(aptot + ZPHptot) , (9) 
Q = q(a PA + 3(3Hp A ) . (10) 



Finally, some brief concluding remarks are given in Sec. |IV 



II. OBSERVATIONAL DATA 



In the present work, we will consider the latest cosmological observations, namely, the 557 Union2 
Type la Supernovae (SNIa) dataset [18], the shift parameter R from the Wilkinson Microwave Anisotropy 
Probe 7-year (WMAP7) data [l9j], and the distance parameter A of the measurement of the baryon 
acoustic oscillation (BAO) peak in the distribution of SDSS luminous red galaxies pfl l2l|. 

The data points of the 557 Union2 SNIa compiled in [l8| are given in terms of the distance modulus 
Pobs(zi). On the other hand, the theoretical distance modulus is defined as 

Pth(zi) = 5 log w D L (zi) + po , (11) 

where p,Q = 42.38 — 51og 10 h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 

D L (z) = {l + z) f ~^— v (12) 
Jo b[z,p) 

in which E = H/Ho, and p denotes the model parameters. Correspondingly, the x 2 from the 557 Union2 
SNIa is given by 

2y _y [Pobs(Zj) - Pth(Zi)} 2 , , 

X ^ (P> ^ a 2 ( Zl ) ' 1 J 

where a is the corresponding \a error. The parameter p,Q is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over /io. However, there is an alternative 
way. Following [2^, , the minimization with respect to po c an be made by expanding the x 2 of Eq. (fl"3"|) 
with respect to po as 



X 2(p) = i-2 Mo B + ^(5, 



(14) 
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where 

i2 



(P) ~ V C>) 



^, x = Aw(z») ~ Hth{zj\y.n = 0,p) (7 = V 



Eq. (|14p has a minimum for //q = B/C at 

= A(p) - ^£ . (15) 

Since xt min = Xu, m i n obviously, we can instead minimize xfi which is independent of (j,q. 

There are some other relevant observational data, such as the observations of cosmic microwave back- 
ground (CMB) anisotropy [l!| and large-scale structure (LSS) [20[. However, using the full data of CMB 
and LSS to perform a global fitting consumes a large amount of computation time and power. As an 
alternative, one can instead use the shift parameter R from the CMB, and the distance parameter A of 
the measurement of the BAO peak in the distribution of SDSS luminous red galaxies. In the literature, 
the shift parameter R and the distance parameter A have been used extensively. It is argued that they 
are model- independent [24j . while R and A contain the main information of the observations of CMB 
and BAO, respectively. 

As is well known, the shift parameter R of the CMB is defined by [24], [25j] 

2 * dz 

m 

where fl m o is the present fractional energy density of pressureless matter; the redshift of recombination 
Z.J, = 1091.3 which has been updated in the Wilkinson Microwave Anisotropy Probe 7- year (WMAP7) 
data [n|. The shift parameter R relates the angular diameter distance to the last scattering surface, the 
comoving size of the sound horizon at and the angular scale of the first acoustic peak in CMB power 
spectrum of temperature fluctuations 0, [25j]. The value of R has been updated to 1.725 ± 0.018 from 
the WMAP7 data pj|. On the other hand, the distance parameter A of the measurement of the BAO 
peak in the distribution of SDSS luminous red galaxies is given by 20] 



R = / , (16) 



1 f Zb dl 



zb Jo E(z) 



(17) 



where z b = 0.35. In 21], the value of A has been determined to be 0.469 (n s /0.98)" - 35 ± 0.017. Here the 
scalar spectral index n s is taken to be 0.963, which has been updated from the WMAP7 data [l!|. So, 
the total x 2 is given by 

X 2 = xl + Xcmb + Xbao > ( 18 ) 

where xl is given in Eq. (J5), Xcmb = ( R ~ Robs) 2 /v 2 R and x 2 B ao = ( A ~ A obs ) 2 /a 2 A . The best-fit 
model parameters are determined by minimizing the total x 2 ■ As in [2(| [22}, the 68.3% confidence 
level is determined by A% 2 = X 2 ~ Xmin — 1-0, 2.3 and 3.53 for n p = 1, 2 and 3, respectively, where 
rip is the number of free model parameters. Similarly, the 95.4% confidence level is determined by 
Ax 2 = X 2 ~ Xmin < 4 A 6 - 17 and 8 - 02 for n p = 1,2 and 3, respectively. 



III. COSMOLOGICAL CONSTRAINTS ON THE SIGN-CHANGEABLE INTERACTIONS 



In this section, we consider the cosmological constraints on the sign-changeable interactions given in 
Eqs. f8j — (fT0| . by using the observational data mentioned in the previous section. 
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FIG. 1: The 68.3% and 95.4% confidence level contours in the Q m o — P plane for the case of Q = 3j3qHp m . The 
best-fit parameters are also indicated by a black solid point. 



A. The case of Q = q(ap m + 3/3Hp m ) 

Firstly, we consider the case of Q — q(ap m + 3/3 Hp m ) given in Eq. ([5]). Substituting it into Eq. ([1]), 
one can find that 

Pro = — - • 3Hp m . (19) 
1 — aq 

Then, substituting into Eq. ([5]), we can finally obtain 

= . 3qHpm . (20) 

1 — aq 



From Eq. ([7]), we have 

Substituting into Eq. (IT^l) . we find that 



p m = ~B. (21) 



H = ^ — - • 3HH , (22) 
1 — aq 

which is in fact a second-order differential equation for H(t). We can change the time t to scale factor a 
with the help of the universal relation / = Haf for any function / (where a prime denotes the derivative 
with respect to scale factor a), and recast Eq. (j22|) as 

aH" + -jtH' 2 + H' = ^— - • 3ff' , (23) 
H 1 — aq 

which is a second-order differential equation for H{a). Note that the deceleration parameter 

q = -l- ^L = -l-- H ' (24) 
is also a function of H and H' . Unfortunately, if a =^ 0, there is no analytical solution for the second-order 
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FIG. 2: Q m , Qa, q and w e g as functions of redshift z with the best-fit parameters for the case of Q = 3j3qHp n 



differential equation (12"B"1) . because one will encounter a transcendental equation. Therefore, we consider 
only the case of a = in this work. In this case, the sign-changeable interaction reads 



Q = 3(3qHp m . 

By solving the second-order differential equation (|2"B1 with a = 0, we find that 

r , „ -i 1/(2+3/3) 

H{a) = Ci 2 3C U (1 + p) - (2 + 3/3) a~^ 1+ ^ 



(25) 



(26) 



where C\\ and C12 are both integral constants, which can be determined in the following. From Eq. PTj) . 
we find that the fractional energy density of dark matter is given by 



a 



K 2 p„ 



3H 2 



Substituting Eq. (gU into Eq. ([27]). we have 

- -7) 



2iY 



2(1 + /3) 



2aiT 
3iJ 



2 + 3/? - 3Cn (1 + /3) a 3 ^) 



Requiring O m (a = 1) = Cl m o, we obtain 



di = 



fi m0 (2 + 3/3)-2(l + /3) 



3O m0 (l + /3) 

On the other hand, requiring H (a = 1) = -ffoj from Eq. (|26[) we can find that 

Cia = Ho [3C U (1 - (2 + 3/3) p 1/(2+3/5) . 
Substituting Eqs. and ([30) into Eq. (gBJ, we finally obtain 



p _H_( 2 + 3/3 

£; -^o-\ 1 -2(rf7)" m0 



l-(l + z) 



3(1+/?) 



1/(2+3/3) 



(27) 



(28) 



(29) 



(30) 



(31) 
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There are two free model parameters, namely fl m o an d f3. Note that when (3 = 0, Eq. (I3T1) reduces to 

E{z) = [O m0 (l + 2) 3 + (1 - O m0 )] V2 , i-e., the one of ACDM model. 

By minimizing the corresponding total \ 2 m Eq. (TT5|) . we find the best- fit parameters J7 TO o = 0.2738 
and (3 = —0.010, whereas Xmin = 542.725. In Fig. [TJ we present the corresponding 68.3% and 95.4% 
confidence level contours in the fl m o — /3 plane. Obviously, the current observational data slightly prefer a 
negative j3. We are also interested to the fractional energy densities il m given in Eq. (j28|) and ft\ = 1— f2 m , 
the deceleration parameter q given in Eq. and the effective equation-of-state parameter (EoS) 

w c ff = ptot/ptot = (2(? — l)/3. We present them as functions of redshift z with the best-fit model 
parameters in Fig. [5] It is easy to find the transition redshift z t = 0.7489 where the universe changes 
from deceleration (q > 0) to acceleration (q < 0). Since the best-fit f3 is negative, dark matter decays 
into dark energy (Q < 0) when z > z tl and dark energy decays into dark matter (Q > 0) when z < z t . 
The interaction Q crosses the non- interacting line {Q = 0) at z t . 



Secondly, we consider the case of Q — q(ap to t + 3f3Hp to t) given in Eq. ([9]). From Eq. ([6]), it is easy to 
find ptot = 3iZ 2 / k 2 . Substituting into Eq. ©, we can finally obtain 



0.14 




0.32 



FIG. 3: The same as in Fig. [T] but for the case of Q = 3f3qHp to t with the condition j3 > 0. 



B. The case of Q — q(ap to t + S/3Hp tot ) 




(32) 



Substituting Eqs. ([2T]) and ([32]) into Eq. ([T]), we have 



H + 3HH (1 + aq) + ^/3qH 3 = . 



(33) 



Similarly, we recast it as 



aH" + —H' 2 + (4 + 3aq) H' + 



9fiqH 
2a 



= 0, 



(34) 
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FIG. 4: The same as in Fig. [2] but for the case of Q = 3f3qHp to t with the condition /3 > 0. 



which is a second-order differential equation for H(a). Note that the deceleration parameter q is also a 
function of H and H' [cf. Eq. (|24l) ]. Similar to the case of Q = q(ap m + 3f3Hp m ), we consider only the 
case of a = in this work. In this case, the sign-changeable interaction reads 



Q = 3f3qHp tot ■ 

By solving the second-order differential equation ([Ml) with a = 0, we find that 

£T(0) = C 2 2 ■ a -3(2-3 / 3+r 1 )/8 . ^3r l/2 + _ 

where C21, C22 are both integral constants, and 



ri = V 4 + P (4 + 9/3) . 
Substituting Eq. ([35J into Eq. (|2"7]). we have 



1 



2-3/3 + 



2C 2 i 



x 3ri / 2 + C-21 



In 



Requiring f2 m (a = 1) = H m o, we obtain 

C 2 i = -1 



2n 



2 - 3/? - 4ft m0 + n ' 
On the other hand, requiring H(a = 1) = Hq, from Eq. (j3"6")l we can find that 

C22 = H Q (i + c 21 y 1/2 . 

From Eqs. (j3"6l and (HtJl) . it is easy to obtain 



£= — = (1 + 2)3(2-3^0/8. 
#0 



(l + z)- 3 ^/ 2 +C7 21 



,1/2 



I + C21 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 
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FIG. 5: The same as in Fig. [TJ but for the case of Q = 3/3qHp to t without the condition /3 > 0. 

where C21 and r\ have been given in Eqs. ([59"]) and (1571) . respectively. There are two free model parameters, 

1/2 

namely tt m o and f3. Note that when (3 — 0, Eq. (|4*Tj) reduces to E(z) = [fi m0 (l + z ) 3 + (1 — ^ ra o)] , 
i.e., the one of ACDM model. 

Imposing the condition < il m < 1 when a — > 0, we have /? > from Eq. (j38|) . Under this condition, 
by minimizing the corresponding total x 2 in Eq- (TT51) . we find the best-fit parameters f2 m o = 0.2701 and 
fi = 0.0, whereas Xmin = 542.919. In Fig. [3J we present the corresponding 68.3% and 95.4% confidence 
level contours in the f2 m o — /? plane. In Fig. U] we also present the f2 m given in Eq. (1551) . f^A = 1 — Q m , 
g given in Eq. ((24]) and w e ff = ptot/ptot = (2<z — l)/3 as functions of redshift z with the best-fit model 
parameters. The universe changes from deceleration (q > 0) to acceleration (q < 0) at the transition 
redshift z t = 0.7549. 

However, the above best-fit model with f3 = is in fact the ACDM model without interaction between 
dark energy and dark matter. So, we would like to give up the condition /3 > 0. This means that in the 
early universe we have fi TO > 1 and then J7a < 0, namely, p\ might be negative. Since the negative energy 
density can arise in quantum field theory (see e.g. [28[ for a good review), it is reasonable to consider 
this possibility. Without the condition f3 > 0, by minimizing the corresponding total x 2 in Eq. (fT5|) . we 
find the best-fit parameters Q m0 = 0.2764 and f3 = -0.0247, whereas xLin = 542.711. In Fig. [5j we 
present the corresponding 68.3% and 95.4% confidence level contours in the f2 m o — j3 plane. Obviously, 
the current observational data slightly prefer a negative (3. In Fig. H2 we also present the fi m given in 
Eq. (|38[) . Q\ = 1 — Q m , q given in Eq. (|24[) and w c g = Ptot/ptot = (2<7 — 1)/3 as functions of redshift z with 
the best-fit model parameters. It is easy to find the transition redshift z t — 0.7688 where the universe 
changes from deceleration (q > 0) to acceleration (q < 0). Since the best-fit (3 is negative, dark matter 
decays into dark energy (Q < 0) when z > z tl and dark energy decays into dark matter (Q > 0) when 
z < z t . The interaction Q crosses the non-interacting line (Q = 0) at z t . 



C. The case of Q = q(apA + 3/3Hpa) 

Finally, we consider the case of Q = q(ap\ + 3(3Hp\) given in Eq. (fTO]) . Substituting it into Eq. ([5]), 
one can find that 

3(3qHp A 

Pa = — 7- ■ I 42 ) 

1 + aq 
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FIG. 6: The same as in Fig. [2| but for the case of Q = 3/3qHp to t without the condition j3 > 0. 

Then, substituting into Eq. (fTO)) . we can finally obtain 

= 3/3gHp A 
1 + aq 

From Eqs. ([5]) and ([7J [or equivalently Eq. (|21l) ]. we have 

3 



k 2 ( 



3# 2 + 2# 



Substituting Eqs. ([21]), 63) and flU} into Eq. ([T]), we find that 

/3q \ 9PqH 3 



H + 3HH 1 



Similarly, we recast it as 



aH" 



H 



-H' 



1 + aqJ 2(1 + aq) 
3(3q \ H , , 9PqH 



= 0. 



1 + aq 



2a(l + ag) 



= 0, 



(43) 



(44) 



(45) 



(46) 



which is a second-order differential equation for H(a). Note that the deceleration parameter q is also 
a function of H and H' [cf. Eq. (|24p]. Unfortunately, if a 0, there is no analytical solution for the 
second-order differential equation (1461) . because one will encounter a transcendental equation. Therefore, 
we consider only the case of a = in this work. In this case, the sign-changeable interaction reads 



Q = SPqHpA . 

By solving the second-order differential equation (|4"6l with a — 0, we find that 

H(a) = C 32 ■ a -3(2-5/ 3+ , 2 )/[4( 2 ^3/ 3 )] . ( a 3r 2 /2 + 

where C31, C32 are both integral constants, and 

r 2 = y/(2-p) 2 = I 2 - f3 \ . 



(47) 



(48) 



(49) 
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FIG. 7: The same as in Fig. \T\ but for the case of Q = 3f3qHpfi 



Substituting Eq. (|48j) into Eq. ([27]), we have 

\l<rn 



2 (2 - 3/3) 

Requiring Cl m (a = 1) = Q, m a, we obtain 

C 3 i = -1 + 



2-5/3 + 



2G 



:S1 



a 3r 2 /2 + Cs] 



-1 r 2 



2r 2 



From Eqs. 



and (I521) . it is easy to obtain 



(1 + z )3(2-5/3+r 2 )/[4(2-3,3)] 



(1 



-3r 2 /2 



1 + C 3 1 



2 - 5/3 + r 2 + 2ft m0 (3/3 - 2) 
On the other hand, requiring H(a = 1) = iJo, from Eq. ([48]) we get 

C 32 =ff (l + C 31 ) 1/(3 ^ 2) . 



C31 



n 1/(2-3/3) 



(50) 



(51) 



(52) 



(53) 



where C31 and r 2 have been given in Eqs. (|51[) and (|49p . respectively. There are two free model parameters, 



Or. 



I 1/2 



namely f2 m o and /3. Note that when /3 = 0, Eq. ((53)) reduces to -E'(z) = [^ m o(l + z ) 3 + (1 
i.e., the one of ACDM model. 

By minimizing the corresponding total \ 2 m Eq. (|18[) . we find the best-fit parameters fl m o = 0.2717 
and /3 = 0.0136, whereas Xmin = 542.778. In Fig. [71 we present the corresponding 68.3% and 95.4% 
confidence level contours in the O m o — P plane. Obviously, the current observational data slightly prefer 
a positive /3. In Fig. [5[ we also present the VL m given in Eq. ([5T))) , = 1 — fl m , q given in Eq. ([M]) 
and !o c (f = ptot/ptot — (2? — l)/3 as functions of redshift z with the best-fit model parameters. It is 
easy to find the transition redshift z t = 0.7398 where the universe changes from deceleration (q > 0) 
to acceleration (q < 0). Since the best-fit /3 is positive, dark energy decays into dark matter {Q > 0) 
when z > Zt, dark matter decays into dark energy (Q < 0) when z < z t . The interaction Q crosses the 
non-interacting line (Q = 0) at z t . 



12 




FIG. 8: The same as in Fig. [2] but for the case of Q = 3/3qHpA- 



IV. CONCLUDING REMARKS 



Recently, Cai and Su [15j found that the sign of interaction Q in the dark sector changed in the 
approximate redshift range of 0.45 <z< 0.9, by using a model-independent method to deal with the ob- 
servational data. In fact, this result raises a remarkable problem, since most of the familiar interactions 
cannot change their signs in the whole cosmic history. Motivated by the work of Cai and Su, we have 
proposed a new type of interaction in a previous work [l6[ . The key ingredient is the deceleration param- 
eter q in the interaction Q, and hence the interaction Q can change its sign when our universe changes 
from deceleration (q > 0) to acceleration (q < 0). In the present work, we consider the cosmological 
constraints on this new type of sign-changeable interactions, by using the latest observational data. We 
find that the cosmological constraints on the model parameters are fairly tight. In particular, the key 
parameter /? has been constrained to a narrow range. 

Some remarks are in order. Firstly, we briefly consider the comparison of these models. For convenience, 
we also consider the well-known ACDM model in addition. In fact, it corresponds to the decaying A model 
with Q = 0. Fitting ACDM model to the observational data considered in the present work, it is easy to 
find the corresponding best-fit parameter £l m0 = 0.2701, whereas Xmin — 542.919. Of course, we would 
like to also consider the decaying A model with a traditional interaction Q = 3f3Hp m which cannot 
change its sign in the whole cosmic history. The corresponding E(z) can be found in e.g. 12|, namely 



E(z) 



fl , ^3(1-/9) , (i _ _^m0_\ 



1/2 



(54) 



Fitting to the same observational data, we find the best-fit parameters f2 m o = 0.2731 and fi = —0.0021, 
whereas x m in = 542.735. A conventional criterion for model comparison in the literature is Xmin/dof, 
in which the degree of freedom dof = N — k, whereas N and k are the number of data points and the 
number of free model parameters, respectively. We present the Xmin/ dof for all the 6 models in Table |TJ 
On the other hand, there are other criterions for model comparison in the literature, such as Bayesian 
Information Criterion (BIC) and Akaike Information Criterion (AIC). The BIC is defined by [29l [3l| 



BIC = -2 In C max + k In N , (55) 
where C max is the maximum likelihood. In the Gaussian cases, Xmin = ~2]n.£ max - So, the difference in 
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Model 


ACDM 


Q = 3f3Hp m 


Q = 3/3qH Pm 


Q = 3/3qHp to t 


Q = 3/3qHptot 


Q = WqHpA 




(Q = o) 






with /3 > 


without /3 > 




Best fit 


a m0 = 0.2701 


n m0 = 0.2731 


n m0 = 0.2738 


« ra0 = 0.2701 


f7 m0 = 0.2764 


Q m0 = 0.2717 






P = -0.0021 


p = -0.010 


/9 = 0.0 


/3 = -0.0247 


/3 = 0.0136 


Xmm 


542.919 


542.735 


542.725 


542.919 


542.711 


542.778 


k 


1 


2 


2 


2 


2 


2 


Xmin /'dof 


0.9730 


0.9744 


0.9744 


0.9747 


0.9743 


0.9745 


ABIC 





6.142 


6.132 


6.326 


6.118 


6.185 


AAIC 





1.816 


1.806 


2.0 


1.792 


1.859 


Rank 


1 


4 


3 


6 


2 


5 



TABLE I: Summarizing all the 6 models considered in this work. 

BIC between two models is given by ABIC = A X 2 mm + AfclniV. The AIC is defined by [13, S3 

AIC = -2ln£ max + 2k. (56) 

The difference in AIC between two models is given by AAIC = Ax^ nin + 2Ak. In TableHJ we also present 
the ABIC and AAIC of all the 6 models considered in this work. Notice that ACDM has been chosen 
to be the fiducial model when we calculate ABIC and AAIC. From Table HI it is easy to see that the 
rank of models is coincident in all the 3 criterions (Xmin/d°f > BIC and AIC). The ACDM model is still 
the best one. However, it is well known that ACDM model is plagued with the cosmological constant 
problem and the coincidence problem (see e.g. [l|). On the other hand, there are some observational 
evidences for the interaction between dark energy and dark matter [l3|, E3| , and the coincidence problem 
can be alleviated in the interacting dark energy models. Therefore, it is still worthwhile to study the 
interacting dark energy models. Although the model with traditional interaction (which cannot change 
its sign) is very close to the other models with sign-changeable interactions, the latter are phenomenally 
richer (see e.g. 16]). Therefore, we consider that the models with sign-changeable interactions deserve 
further investigations. 

Secondly, we note that the case of Q = 3f3qHpA is fairly different from the cases of Q = 3f3qHp m and 
Q = 3f3qHp to t- Comparing Fig. [7] with Figs. [TJ|3] and [5l it is easy to see that the direction of contours in 
the fi m o — /3 plane is rightward for the case of Q = 3(3qHp\, whereas the ones are leftward for the cases 
of Q = 3/3qHp m and Q = 3[3qHp tot . From Table HI we find that the best- fit j3 is positive for the case of 
Q = 3/3qHp\, whereas the ones are negative (or zero) for the cases of Q — 3f3qHp m and Q — 3f3qHp to t- 
This means that in the case of Q = 3(3qHpA the interaction Q crosses the non-interacting line (Q = 0) 
from above to below, whereas in the cases of Q = 3(3qHp m and Q = 3(3qHp to t the interaction Q crosses 
the non-interacting line (Q = 0) from below to above. This is physically interesting, because Q > 
means that the energy transfers from dark energy to dark matter, whereas Q < means that the energy 
transfers from dark matter to dark energy. 

Finally, in this work the role of dark energy is only played by the decaying A (vacuum energy) , whereas 
the parameter a in the sign-changeable interactions are chosen to be zero. So, the constraints obtained 
in this work cannot be directly used to the models different from the ones considered here. In fact, the 
interacting dark energy models with sign-changeable interactions can be generalized. For instance, one 
can choose dark energy to be the one with a constant or variable EoS (including parameterized EoS, 
or even the ones of quintessence, phantom, k-essence, Chaplygin gas, quintom, hessence, holographic 
or agegraphic dark energy, and so on). Of course, one can also let the parameter a be free and then 
constrain the models numerically. 
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